Metastability and the Casimir Effect in Micromechanical Systems 
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Electrostatic and Casimir interactions limit the range of positional stability of electrostatically- 
actuated or capacitively-coupled mechanical devices. We investigate this range experimentally for 
a generic system consisting of a doubly-clamped Au suspended beam, capacitively-coupled to an 
adjacent stationary electrode. The mechanical properties of the beam, both in the linear and 
nonlinear regimes, are monitored as the attractive forces are increased to the point of instability. 
There "pull-in" occurs, resulting in permanent adhesion between the electrodes. We investigate, 
experimentally and theoretically, the position-dependent lifetimes of the free state (existing prior to 
pull-in). We find that the data cannot be accounted for by simple theory; the discrepancy may be 
reflective of internal structural instabilities within the metal electrodes. 
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The technology of micro electro mechanical systems 
(MEMS) now routinely allows fabrication of miniature 
movable structures on-chip. This new branch of mi- 
croelectronics enables the realization of a variety of 
miniature, fully integrated sensors and actuators, with 
a rapidly growing range of applications. One common 
building block within MEMS is the capacitively-tunable 
resonator (CTR), comprising a movable electrode that 
serves as a mechanical resonator, and an adjacent sta- 
tionary electrode (see jlj and references therein). The 
capacitive coupling between the electrodes allows electro- 
static control of both the mean position of the resonator 
and its resonance frequency. In addition to the elec- 
trostatic interaction, such electrodes are coupled via the 
Casimir effect ||. This quantum electrodynamical 
effect originates from the dependence of the ground state 
energy of the electromagnetic field upon boundary con- 
ditions and leads to an observable forces between macro- 
scopic bodies. Due to its relatively short range, this force 
has only a small effect on the dynamics of macroscopic 
mechanical systems. However, the Casimir force can 
play a major role in modern MEMS where typical dis- 
tances between neighboring surfaces can be on the sub- 
micron length scale These attractive electrostatic 
and Casimir forces give rise to a mechanical metastabil- 
ity that may cause stiction (for a review see Q). In 
the present case it is manifested as collapse (pull-in) of 
the movable resonator onto the nearby electrode, result- 
ing in their permanent adhesion. This phenomena can 
be a principal cause of malfunctioning in MEMS. More- 
over, it limits the range of stable operation of CTRs. 

In the present work we study experimentally and theo- 
retically the range of tunability and mechanical metasta- 
bility in a CTR system made of a doubly clamped Au 
beam separated from an adjacent counterelectrode by 
a small vacuum gap. Once the beam is brought to 
contact with the electrode, permanent adhesion occurs. 
This indicates that the free state of this system is merely 
metastable, and that the state of contact (after pull-in) 
has lower energy due to the strongly attractive Casimir 



interaction. The potential barrier separating these two 
states determines the lifetime of the free state. The 
barrier can be reduced by introducing electrostatic at- 
traction, namely by applying a DC voltage, V dc between 
the beam and the electrode. In our measurements we 
study the mechanical properties of the beam as Vdc is 
gradually increased to the point of pull-in. Near this 
critical point all of the measured properties of the beam 
show strong dependence on Vdc- Using a simple model 
we calculate the shape of the potential barrier confin- 
ing the metastable mechanical state, and then estimate 
the escape rate via both thermal excitation and quantum 
tunneling. 

Fig. 1(a) is a micrograph displaying a side view of 
our device. The details of the fabrication process are 
given elsewhere 0. The structure is designed to al- 
low full characterization of the beam's properties to yield 
straightforward and unambiguous interpretation of our 
results. We use bulk (rather than surface) micromachin- 
ing, which allows the substrate to be completely removed 
beneath the sample. This greatly simplifies the bound- 
ary conditions of the electromagnetic field in the vicinity 
of the sample. Moreover, we avoid using multilayered 
structures, since their internal stresses are generally im- 
portant, and are difficult to model theoretically. The 
beam has length I = 200 /im, width a = 0.28/xm and 
thickness t = 0.25/im (measured using an atomic force 
microscope). At its center, the stationary electrode has 
a 20/im long rectangular protrusion. This is separated 
from the beam by a vacuum gap of g = 1.28 /im. 

All measurements are done at room temperature, in- 
situ within a commercial scanning electron microscope 
(SEM). A voltage V = V dc + V ac cos (2nvt) is applied be- 
tween the beam and the electrode (see Fig. 1(b)). The 
static response of the beam is measured by operating the 
SEM in scanning mode, imaging the device, and using 
digital image processing to extract experimental data. 
The dynamical response is measured with a stationary 
electron beam focused at a point near the edge of the 
Au beam. To detect mechanical displacement we use a 
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spectrum analyzer to monitor the output signal from the 
photomultiplier serving as a secondary electron detector. 
Note that this detection scheme is sensitive almost ex- 
clusively to motion in the plane of the sample. 

To characterize the mechanical properties of the beam, 
namely stress and stiffness, we measure its resonance fre- 
quencies. For small stiffness they are given by: 



v n = nv Q [1 + 2C + (4 + nV/2) ( 



(1) 



where C 2 = Ea 3 t/12Tl 2 , v = y/T/pA/2l, with E be- 
ing Young's modulus, T is the stress, and p is the mass 
density , || . The dimensionless parameter £ indicates 
the relative effect of stiffness compared to stress on the 
dynamics of the beam. Note that the terms that make 
the spectrum unequally spaced are of order O (£ 2 ) . 

We excite the beam by applying AC voltage (Vd c = 0) 
and find three modes with frequencies vi = 185.53 kHz, 
v 2 = 372.4 kHz, and v z = 563.8 kHz. Using Eq. (g) 
we find £=0.017±0.002, thus stiffness has only a small 
effect on the dynamics of the beam. In what follows, we 
ignore altogether corrections due to stiffness. The stress 
is found from the measured resonance frequencies to be 
T = 7.4 x I(T 6 N. 

Next we gradually increase Vdc and measure the me- 
chanical properties of the beam until pull-in occurs. Fig- 
ure 2(a) shows the normalized static deflection of the 
beam's center, X\/g, as a function of Vdc- After com- 
pleting the measurements at Vdc = 29 V, pull-in occurred 
with this particular beam. Simultaneously applying the 
small AC excitation with the DC bias allows us to mea- 
sure the frequency of the fundamental mode v\ = u)i/2n. 
Figure. 2(b) shows the ratio — Aw/wo = — (^i — ojq) /u>q, 
where ujq is the angular frequency at Vdc = 0. Both X\ 
and uj\ show strong dependence on Vdc near pull-in. 

Further insight can be gained by studying the dynam- 
ics of the beam in the nonlinear regime. For oscilla- 
tions with a large amplitude the harmonic approximation 
breaks down and effects due to higher order terms of the 
potential are observable. Figure 3 shows the measured 
response of a beam as a function of frequency around the 
fundamental resonance for a discrete range of ac drive 
{Vac) levels, at fixed DC bias, Vdc = 20 V. As V ac is in- 
creased, the response becomes asymmetric, with a sharp 
drop occuring on the high frequency side. Moreover, a 
strong hysteresis is found for high V ac when the frequency 
is scanned back and forth across the resonance (see inset 
of Fig. 3). Theoretically, the response function in the 
nonlinear regime is given to lowest order by || : 
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where / is the amplitude of the external force, a is the 
amplitude of oscillations, m is the mass, and k is a param- 
eter describing the amplitude-dependent frequency shift 
due to nonlinearity. 



From the shape of the resonance peaks in the nonlinear 
regime we extract the nonlinear frequency-shift parame- 
ter, k, as a function of Vdc (Fig. 2(c)). For relatively 
small Vdc we find k > 0. For higher Vdc K becomes nega- 
tive, and close to the point of pull-in k drops quite rapidly 
as a function of Vdc, indicating the strong nonlinearity of 
the potential in this regime. 

We describe our system using a simple one dimensional 
model as proposed in Ref. (lO) (see Fig. 4(a)). This 
model employs rather simplistic approximations and thus 
cannot be expected to provide accurate quantitative pre- 
dictions - it is extremely useful, however, in providing 
physical insight and a qualitative description of the dy- 
namics. A mass m coupled to a string with a constant 
k = muj 2 represents the doubly-clamped beam. An ad- 
jacent surface at distance g interacts with the mass elec- 
trostatically, as well as by the Casimir effect. For our 
device geometry it is appropriate to model these forces, 
F, as applied locally to the beam's center. In this case 
the displacement of the beam's center is Fl/AT, and the 
effective spring constant is thus k — AT /I. 

The total potential energy is given by: 

U t (x) = U E (x) + U v (g-x) + U c (g ~ x) . (3) 

The first term Ue (s) = (1/2) ks 2 is the beam's elastic po- 
tential. The second term is the electrostatic potential in- 
troduced by applying a DC voltage, V, between the beam 
and its counterelectrode. We estimate the capacitance 
C between the beam and the electrode using a simple 
parallel plate approximation, thus Uy (s) = —e^AV 2 /2s, 
where eo is the vacuum permittivity and A is the area. 
The third term is the Casimir potential. As in the elec- 
trostatic case, we employ a parallel plate approximation 
leading to U c (s) = (-7r 2 /720) (hcA/s 3 ), where h is the 
Planck's constant and c is the velocity of light || . Using 
these expressions we rewrite (/ ( in a dimensionless form: 



U t = U E (g) 



ft 



1 - x/g 1 - ( X /gf 



(4) 



where a = —Uv{g) /Usig) — eoAV 2 /kg 3 represents 
the relative importance of the electrostatic interaction 
(compared with the elastic term), and similarly (3 = 
-U c {g) /U E (g) = (tt 2 /360) (hcA/kg 5 ) for the ratio of 
Casimir interaction to the elastic term. 

For relatively small values of a and ft (see figure 4(b)) 
the potential energy has a local minimum at x = x\ 
associated with a metastable state. Figure 4(c) shows 
a typical example of Ut in the metastable regime. The 
angular frequency of small oscillations around the local 
minimum is u\. The barrier that separates the local 
minimum at x = X\ from the global minimum at x = 
g has height Ub- Below, we calculate the parameters 
toi and Ub for two special cases relevant to the present 
experiment. 
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Consider first the case of a « 1 and <C 1. In 
this case we find Xx/g — (a + 3/3)/2 and lji/uq = 
(1 — a/2 — 3/3). Since the escape rate for this case is 
very small, we do not estimate the associated Ub- 

In the next case we consider the instability occurring 
with a strong electrostatic and a weak Casimir interac- 
tion. Note first that if = then a local minimum and a 
potential barrier exist only for a < 8/27. At the critical 
point, a = 8/27, the local minimum and the local maxi- 
mum points coincide at x/g = 1/3. We therefore expand 
the potential near that point, assuming e = 8/27 — a <C 1 
(e > 0) and also <C 1. To lowest order we find 



xx/g = 1/3 - (1/8) ^32e - 216/3 
ux/w Q = (3/4) (32e - 216/3) 1/4 
U b /kg 2 = (3/512) (32e - 216/3) 3/2 



(5a) 
(5b) 
(5c) 



The critical point is thus shifted to a = 8/27 — 27/3/4 due 
to the Casimir effect. The nonlinear frequency-shift pa- 
rameter k is calculated from the coefficients of the Taylor 
expansion of the potential around the point x = xx |§- 
To lowest order we find ng 2 oc (32e — 216/3)~ 3 ^ 4 . 

To compare the above results with theory we substi- 
tute the area, A, with an effective area, A e ff = 2.3A, 
chosen to ensure that the measured and calculated criti- 
cal points coincide. This substitution serves to account, 
in part, for the crude approximations employed in the 
model. The solid lines in Fig. 2(a) and 2(b) show the 
theoretical results of Eq. ((5a]). Agreement with experi- 
ment for both xx and u>x is moderately good despite the 
simplicity of the model. Note that our measurements 
within the nonlinear regime are strongly affected by non- 
linearity in the detector's response. Given that this is 
not properly accounted for in our model, we do not at- 
tempt comparison with theory for the parameter k. 

We now turn our attention to the lifetime of the 
metastable mechanical state before pull-in. At finite 
temperatures, thermal excitation above the energy bar- 
rier can allow the system to "escape" from the freely- 
moving to the pulled-in state. The rate for this generic 
type of process was first found by Kramers |l2) , via the 
Fokker-Planck equation. The thermal escape rate for 
the case of small damping is 



Tt= Qk^f exp 



Ub \ 
k B Tj 



(6) 



This result is valid when U b /k B T > 1 and U b /k B T > 
1/Q. Alternatively, at low temperatures the system es- 
capes via quantum tunneling |]ll"| . Caldeira and Leggett 
(CL) H H have explored what is termed macroscopic 
quantum tunneling (MQT ) in systems with many degrees 
of freedom. In the case of small damping they find that 
the MQT rate is 
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Here £i , £2 and £3 are positive constants of order unity. 
The rate for the case of no dissipation (1/Q = 0) co- 
incides with the conventional WKB result for tunneling 
within a noninteracting, single-particle picture. The CL 
model yields a description of the effect of dissipation upon 
MQT involving only a single parameter, Q, that sub- 
sumes all of the precise details of the interaction between 
the beam and the environment. Equation (^) explicitly 
demonstrates that the effect of dissipation is to suppress 
the tunneling rate. 

In the present experiment the temperature is relatively 
high (k B T 3> huti) therefore the rate of thermal escape 
is much larger compared to the rate of quantum tunnel- 
ing. However, the opposite situation, namely Tq > T t , 
may be realized with present-day technology. For high 
frequency mechanical oscillators (see for example [15)) 
the quantum limit, k B T < Hloi, is reachable at cryogenic 
temperatures. The model employed in the present pa- 
per may be useful in designing an experiment to study 
phenomena in this regime. 

Comparison between the simple theory and experi- 
ment indicates significant disagreement. We extract 
the parameter 32e — 216/3 from Eq. (5b) and the value 
u>x/wo — 0-81 found experimentally when the pull-in oc- 
curs. This parameter is substituted in Eq. (5c) to find 
Ub- This yields Ub/k B T ~ 10 7 , a value which, when 
substitited into Eq. (^), yields a lifetime that is essen- 
tially infinite for all practical purposes. Experimentally, 
however, we observed pull-in after only a few minutes, 
in striking contradiction with this theoretical estimate. 
Similar devices, having g =1.29, 0.71 and 0.65 /mi, exhib- 
ited pull-in when oji/wo = 0.84, 0.987, and 0.995 respec- 
tively. In these cases, again, the simple theory is far from 
explaining what is experimentally observed. We note, in 
general, that any theory predicting Ub = Akg 2 (uji/ujq) B 
near the critical point, where both A and B are of order 
unity, will exhibit similar discrepancy. This is due to the 
exponential dependence in Eq. ([|), in concert with the 
fact that kg 2 /k B T ~ 10 7 > 1. 

Our results do not appear to be anomalous. We have 
analyzed data from previous experiments involving CTR 
devices made from stress-free silicon 0, in which pull- 
in was also observed upon application of a DC voltage 
between electrodes. Calculating the expected lifetime 
for the experimental conditions of Ref. [[TJ using Eq. (^), 
we find lifetimes that are in similarly gross contradic- 
tion with experimental observations. In principle, slow 
processes of stress relaxation in both cases may be the 
underlying mechanism inducing early pull-in fig] . Fur- 
ther investigations are warranted. In our work we have 
taken care to rule out spurious effects from environmen- 
tal electrical noise sources that could, in principle, drive 
the system out of the metastable state. Measured elec- 
trical noise was far smaller than the steady bias voltages 
applied. Also, "early" pull-in occurred even when the 
electron beam was blanked, ruling out the possibility that 
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some form of parasitic local charging may be operative. 

The present work shows that MEMS can provide an 
ideal approach to studies of mechanical metastability. 
Further experimental and theoretical work will elucidate 
the underlying mechanisms governing the instability of 
these systems. Such work is important not only from a 
theoretical point of view, but also from the standpoint 
of practical device engineering. Moreover, further en- 
hancement of experimental techniques will enable new 
studies of macroscopic quantum tunneling in mechanical 
systems. 
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FIG. 1. (a) Side view micrograph of the device, (b) Exter- 
nal force is applied by a voltage bias and motion is detected 
by focusing the e-beam on a point on the edge of the structure 
and measuring the secondary reflected electrons. 

FIG. 2. (a) Displacement of the center of the beam xi di- 
vided by g. (b) Frequency shift Auu divided by the unper- 
turbed frequency loo- (c) Nonlinear parameter Kg 2 . 

FIG. 3. The mechanical response as a function of frequency 
around the fundamental resonance with Vdc — 20 V and 
Vac =25, 50, 75,.. .,225 mV. The inset shows the hysteresis 
loop for Vac = 225 mV when the frequency is scanned forth 
and back across the resonance. 

FIG. 4. (a) A one dimensional model describing our sys- 
tem, (b) the metastable regime in the a/3 plane, (c) The 
potential energy in the metastable regime has a barrier Ub 
separating the metastable state at x — x\ and the contact 
state at x — g. 
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